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. A model of electromagnetic interaction of relativistic heavy ions 

CN \ which makes it possible to calculate the inclusive and exclusive char- 

acteristics of such interactions is presented. The multistep model in- 
Q\ . eludes (1) the single and double virtual photon absorption by a nu- 

cleus, (2) intranuclear cascades of produced hadrons and (3) statisti- 
cal decay of excited residual nuclei described by evaporation, fission 
and multifragmentation mechanisms. Components of the model were 
^ I verified by comparison of calculation results with experimental data 

obtained with monoenergetic photons and heavy ions at KEK-Tanashi 
Synchrotron and CERN SPS accelerator, respectively. 

X 
b . 

1 Description of the model 

The relativistic Coulomb excitation of nuclei is a well-known phenomenon Q . 
It is expected to play an important role in collisions of ultrarelativistic heavy 
nuclei at new colliders, RHIC and LHC 0, |J. The process may be described 
in the framework of the Weizsacker- Williams method of equivalent photons. 
According to this method the impact of Lorentz-boosted Coulomb field of a 
nucleus on the collision partner may be treated as the absorption of equiva- 
lent photons. Due to the coherent nature of photon emission by the nucleus, 
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such photons are almost real, Q 2 < 1/R 2 , R is the nuclear radius. Therefore 
photonuclear reaction data obtained in experiments with real monoenergetic 
photons and appropriate theoretical models describing photonuclear reac- 
tions may be used to predict the properties of electromagnetic dissociation 
process. 

Since an ultrarelativistic nucleus moving at the impact parameter b with 
the velocity /3 = t>/c~l, (7 3> 1) spends a short time At near the collision 
partner, the virtual photon spectrum contains all the frequencies up to the 
maximum energy E™ ax ~ I /At ~ 7/-R: 

N(E„b) = ^^(^ 2 (x) + ^ 2 (x)), (1) 

where a is the fine structure constant, K and K\ are the modified Bessel 
functions and x = E^b/^fi). The mean number of photons absorbed by the 
collision partner of mass A p is defined by: 

j^max 

m(b)= J N(E 7 ,b)a Ap (E 7 )dE 7 , (2) 

f^mirt 

7 

where a Ap (E 7 ) is the appropriate total photoabsorption cross section ob- 
tained starting from the photoneutron threshold at E™ m ~ 7 MeV. 

Assuming the Poison distribution for the multiphoton absorption with the 
mean multiplicity m{b), Eq.(|2]), the integral cross sections are calculated for 
the first- and second-order processes @] for a particular dissociation channel v. 

j^max 

^ = J dE 1 N^{E 1 )a Ap {E 1 )fP{E 1 ), (3) 

fpmin 
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J J dE l dE 2 N^\E ll E 2 )a Ap (E l )a Ap (E 2 )fi 2 \E ll E 2 ), (4) 



TPmin TPmin 
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with the following single and double photon spectral functions: 

00 

N^{Ex) = 2tt / bdbe- m{b) N{E u b), 



OO 



N {2 \E U E 2 ) = 71 J bdbe- m ^N(E 1 ,b)N(E 2 ,b). 
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The values (Ei) and f\ (Ei,E 2 ) denned as the branching ratios for the 
considered disintegration channel i in the single and double photon absorp- 
tion, respectively, have to be calculated by our multistep model. Similar 
expressions may be written for the differential distributions of produced par- 
ticles on the rapidity, transverse momentum and other variables. 

In order to obtain and several mechanisms are included in a new 
Relativistic ELectromagnetic Dissociation (RELDIS) code for the Monte 
Carlo simulation, namely, the intranuclear cascade of fast particles pro- 
duced after the photon absorption on a nucleon or nuclear pair ||, and the 
evaporation of nucleons and lightest fragments, binary fission or multifrag- 
mentation || at a later stage of interaction. The multifragmentation pro- 
cess dominates when the excitation energy of residual nucleus, E*, exceeds 
3 — 4 MeV/nucleon and is described by the statistical multifragmentation 
model (SMM) ||. For fissile nuclei the evaporation may take place before or 
after the fission. The competition of evaporation and fission is also described 
with the SMM. The lightest fragments may be also created via coalescence 
of fast nucleons into d, t, 3 He or 4 He 0. The decay of highly excited light 
residual nuclei with A < 16 is treated by the Fermi break-up mechanism ||. 
Other details of the calculation scheme may be found in Refs. || [| 

2 Emission of nucleons and lightest fragments 

Depending on the photon energy, E^, and mass number, A p , different pro- 
cesses take place in the nuclear photoabsorption. Due to the contributions 
of several mechanisms: 7iV — > ttN, 'yN — > 2nN and j(np) — > np the 
calculated double differential spectra d 2 a/dQdP of pions and protons have 
complex shapes above the pion production threshold || . The spectra of fast 
particles 7r + ,7r _ ,^ and p predicted by the model were compared in Ref. || 
with available sets of experimental data at 140 < E 1 < 1000 MeV including 
the data obtained with KEK-Tanashi, the 1.3-GeV Electron Synchrotron. A 
satisfactory description of the spectra was obtained. 

Fast hadrons produced after the photon absorption initiate a cascade of 
subsequent collisions with the intranuclear nucleons leading to the heating 
of a residual nucleus. Later the nucleus undergoes de-excitation by means of 
the emission of nucleons and fragments. Because of a low Coulomb barrier 
in light nuclei the rates of proton and neutron emission are comparable. 

Fragment spectra in the photoabsorption on a carbon nucleus are given 
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in Fig. |I|. The low energy part of the deuteron spectra is explained by the 
explosive Fermi break-up while the high energy part is attributed to the 
coalescence mechanism. Since the main part of fast nucleons is emitted in 
the forward direction, the coalescence contribution dominates at small angles. 
On the contrary, the distribution of deuterons from Fermi break-up is nearly 
isotropic. Because of the option to accelerate also light oxygen ions at RHIC, 




Figure 1: Deuteron emission in photoabsorption on carbon. Points: KEK- 
Tanashi data || . Solid histograms: results of the intranuclear cascade calcu- 
lations with coalescence and Fermi break-up. Contributions from coalescence 
mechanism are shown by the dashed, dotted and dash-dotted histograms for 
different values of coalescence parameter [[F]: p = 90,129 and 200 MeV/c, 
respectively. 

the calculations of the electromagnetic dissociation contribution for such ions 
should take into account both of the considered mechanisms. 

The photoabsorption scenario is very different for heavy nuclei, gold and 
lead. In this case neutrons are the most abundant particles produced in 
the electromagnetic collisions of ultrarelativistic nuclei Q. When a heavy 
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nucleus absorbs photons in the Giant Resonance region, 6 < E n < 30 MeV, 
the evaporation model may be used with the assumption E* = E 7 which 
leads to the emission of one or two neutrons. 

On the other side of the equivalent photon spectrum, when E^ reaches 
the value of several GeV, the multiple pion photoproduction on intranuclear 
nucleons becomes the main absorption mechanism. In this case up to 95% 
of the photon energy is released in the fast particles on average. Never- 
theless, the remaining energy deposited in the residual nucleus is sufficient 
for evaporating many neutrons or even multifragmentation ||. As it was 
found in Ref. g for PbPb collisions at SPS (1584 GeV beam) and LHC 
(2. 754 + 2.754 TeV beams) the mean neutron multiplicities are 4.2 and 8.8, 
respectively. The same value for AuAu collisions at RHIC (1004 + 1004 
GeV) is equal to 7.2. 

The neutron multiplicity distributions are shown in Fig. £| They are 
strongly peaked at the single neutron emission channel due to GDR decay, 
while there is a long tail of the multiple neutron emission originating from 
high excitations. These results provide important information for designing 
large-rapidity detectors and zero-degree calorimeters at RHIC and LHC. 




Figure 2: RELDIS predictions for multiplicity distributions of neutrons in 
the electromagnetic dissociation of Pb nuclei at LHC and SPS energies (solid 
and dotted lines, respectively) and Au nuclei at RHIC energies (dashed lines). 
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3 Charge changing reactions at SPS 



CERN SPS accelerator currently supplies the highest energy available for 
heavy ion studies. A detailed understanding of fragmentation mechanisms 
at SPS energies provides an important point for the extrapolation of the 
theoretical and experimental results to RHIC and LHC energies. 




Figure 3: RELDIS results for fragmen- 
tation charge changing cross sections 
of 208 Pb ions at SPS. Points: 
mental data 



experi- 



1C 



Figure 4: RELDIS predictions for 
mass distributions of isotopes with 
charges Z = 75, 76, ...,82 produced 
in the electromagnetic dissociation 
of 208 Pb ions at SPS. 



Partial break-up charge changing cross sections (—7 < AZ < 1) for 158 
AGeV 208 Pb ions interacting with various targets, from H to Pb, were mea- 
sured recently [fTofl. Although the contributions to fragmentation due to the 



electromagnetic and nuclear forces were not separated in the experiment [10 
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one can use RELDIS code to evaluate the role of the electromagnetic dissoci- 
ation for each particular dissociation mode. As it is shown in Fig. |3|, the main 
part of charge changing cross sections in PbPb collisions may be explained 
by the electromagnetic interaction. 

Calculations show that multiple neutron emission also plays a role in 
each of the partial charge changing channels with AZ = —6, — 1 for PbPb 
collisions, see Fig. [|. The lost of each proton is accompanied with a high 
probability by the multiple neutron emission. For example, up to 20 neutrons 
may be emitted in AZ = —2 fragmentation channel, and the emission of 8-10 
neutrons is the most probable process. 

This is evident from the fact that due to a high Coulomb barrier of Pb 
nucleus, photoemission of protons takes place well above the GDR region, 
particularly in the region of quasideuteron absorption, where the excitation 
energy, E* ~ 20 — 50 MeV, is sufficient to evaporate many neutrons. On the 
contrary, the mass distribution for AZ = channel is completely different 
since in this case In and In emission in GDR region strongly dominates. 
However, a long tail of the multiple neutron emission also exists in this 
channel. 

Our expectation of the intense neutron emission in the charge changing 
electromagnetic collisions of heavy nuclei is very different from a simple pic- 



ture assumed in Ref. [IT], where the rates of neutron and proton emission 
above the photoproton threshold were determined by N/ Z ratio of the nucleus 
which undergoes fragmentation. A detailed understanding of the fragmenta- 
tion mechanism may be obtained by measurements of fragment masses and 
detection of neutrons produced in collisions of ultrarelativistic heavy ions. 



4 Conclusions 

A variety of fragmentation mechanisms takes place in the electromagnetic 
dissociation of ultrarelativistic heavy ions, namely the coalescence, Fermi 
break-up, evaporation, fission and multifragmentation. A partial or complete 
disintegration of the colliding nuclei is possible despite the absence of geomet- 
rical overlap of the nuclear densities. Comparison with CERN SPS data on 
PbPb collisions demonstrates the dominance of electromagnetic dissociation 
in the partial charge-changing fragmentation channels with AZ = —6, — 1. 
RELDIS model predicts also a very intense neutron emission in such chan- 
nels. 
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